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ABSTRACT 

We report results from a comprehensive follow-up observing campaign of the 
afterglow of GRB 030226, including VLT spectroscopy, VLT polarimetry, and 
Chandra X-ray observations. In addition, we present BOOTES-1 wide-field ob- 
servations at the time of the occurrence of the burst. First observations at ESO 
started 0.2 days after the event when the GRB afterglow was at a magnitude of 
R ~19 and continued until the afterglow had faded below the detection thresh- 
old (R > 26). No underlying host galaxy was found. The optical light curve 
shows a break around 0.8 days after the burst, which is achromatic within the 
observational errors, supporting the view that it was due to a jetted explosion. 
Close to the break time the degree of linear polarization of the afterglow light was 
less than 1.1%, which favors a uniform jet model rather than a structured one. 
VLT spectra show two absorption line systems at redshifts z = 1.962 ±0.001 and 
at z — 1.986 ± 0.001, placing the lower limit for the redshift of the GRB close 
to 2. We emphasize that the kinematics and the composition of the absorbing 
clouds responsible for these line systems is very similar to those observed in the 
afterglow of GRB 021004. This corroborates the picture in which at least some 
GRBs are physically related to the explosion of a Wolf-Rayet star. 

Subject headings: Gamma rays: bursts - Stars: Wolf-Rayet 



1. Introduction 

Understanding the nature of GRB progenitors remains a primary focus in GRB research. 
The spectroscopic identification of supernova (SN) light underlying the afterglows of the 
nearby bursts 030329 (Hjorth et al. 2003a; Kawabata et al. 2003; Matheson et al. 2003; 
Stanek et al. 2003) and 031203 (Malesani et al. 2004) convincingly demonstrated that core- 
collaps supernovae are physically related to long-duration GRBs. In most cases, however, 
GRB-SNe are too faint at their suspected maximum light (cf. Zeh et al. 2004) to confirm 
their appearance via spectroscopic observations. Revealing the nature of a GRB progenitor 
by indirect methods thus remains an important approach, in addition to an understanding of 
the afterglows themselves. In this context of particular interest is the fine structure detected 



1 Based on observations collected at the European Southern Observatory, La Silla and Paranal, Chile 
(ESO Programmes 70.D-0523 and 70.D-0531). 
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in two of the best monitored afterglow light curves (GRB 021004: e.g., Bersier et al. 2003; 
GRB 030329: e.g., Lipkin et al. 2004) and the kinematics of absorbing clouds revealed in 
early-time spectra of the afterglows of GRB 020813 and 021004 (for a discussion see Chevalier 
et al. 2004; Mirabal et al. 2003; and Schaefer et al. 2003), since both features could be 
signatures of the physical conditions in the GRB environment. Here we report the results of 
our comprehensive observing campaign of the afterglow of GRB 030226. Even though the 
afterglow light curve started faint and faded rapidly, which did not allow us to monitor its 
evolution in great detail, our early-time spectra do reveal strong spectral similarities between 
the afterglow of GRB 030226 and those of GRB 020813 and 021004. 

GRB 030226 was discovered by the HETE-2 satellite on 2003 February 26.15731 UT 
(3:46 UT) with a duration in excess of 100 seconds and a peak flux and fluence in the 30-400 
keV band of ~1.2xl0~ 7 erg cm" 2 s _1 and 5.7xl0~ 6 erg cm -2 , respectively (Suzuki et al. 
2003). The optical counterpart was detected ~2.5 hours later (Fox et al. 2003, 2004) at a 
magnitude R ~ 18.5 (Garnavich et al. 2003) at coordinates (J2000) R.A. = ll h 33 m 04!92; 
Decl. = +25° 53' 55"6 (Price et al. 2003a; von Braun et al. 2003), well within the HETE-2 
Wide Field X-ray monitor (WXM) and Soft X-ray Camera (SXC) error boxes. No optical 
emission coincident in time with the GRB prompt event was detected (Castro-Tirado et al. 
2003; § 4.2). A subsequent Chandra observation led to the detection of the X-ray afterglow 
at a position consistent with that of the optical transient, thus confirming its association 
with GRB 030226 (Pedersen et al. 2004a). 



2. Observations and data reduction 

Imaging of the central part of the SXC error circle in J s , H, K s began ~4.5 hrs after the 
burst using VLT/ISAAC at ESO Paranal. Simultaneously, imaging in the optical in BVR 
was performed using the ESO NTT telescope at La Silla equipped with the multi-purpose 
instrument EMMI. Also, the 3.8-m United Kingdom Infra-Red Telescope (UKIRT) at Mauna 
Kea, Hawaii, started observing the GRB afterglow 6 hours after the burst using the UKIRT 
Fast-Track-Imager (UFTI). Further imaging was performed during the following nights with 
the ESO NTT telescope equipped with the SUperb Seeing Imager (SuSI2). Additional 
data were obtained with the Wendelstein 80-cm telescope near Munich, Germany (using the 
direct imaging CCD camera MONICA), the 2.5-m Nordic Optical Telescope (NOT) equipped 
with the Andalucia Faint Object Spectrograph and Camera ALFOSC, the 1.0-m Jacobus 
Kapteyn Telescope (JKT), and the Italian 3.5-m Telescopio Nazionale Galileo TNG, using 
the focal reducer instrument DOLORES (all at La Palma, Spain), the Tautenburg 1.34- 
m Schmidt telescope equipped with prime focus 4k x 4k CCD camera, the Asiago 1.8-m 
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telescope equipped with the AFOSC camera, UKIRT, and the 8-m Gemini South telescope 
at Ccrro Pachon, Chile. The data base, and the derived magnitudes (not corrected for 
Galactic extinction; E(B — V) ~ 0.02 mag, Schlegel et al. 1998) are summarized in Table 2. 
These magnitudes supersede those obtained by this collaboration and published in the GCN 
Circular archive. 

Regular observations of the sky were also acquired by the wide-field cameras of the 
BOOTES project (Castro-Tirado et al. 1999) at the BOOTES-1 astronomical station in 
Mazagon, Spain. The GRB 030226 location was serepinditously imaged in a simultaneous, 
unfiltered frame (180 s exposure time) taken on 03:45:16 - 03:48:16 UT. Additional images 
were available prior to and following the event. 

2.1. Photometry 

After the normal image processing steps (flatfielding, stacking), all stars were extracted 
using DAOPHOT as implemented in IRAF 2 . Point-spread-function (PSF) fitting was used 
to derive the brightnesses of the afterglow and comparison stars. PSF stars were carefully 
selected by first inspecting the deeper, higher resolution images to ensure that no background 
contamination was present and that all objects were stellar. After extraction, ensemble dif- 
ferential photometry was performed using comparison stars chosen from the field photometry 
file provided by Henden (2003). Typically four to six comparison stars were used, depending 
on the depth of the image and the bandpass. No transformations were made to the standard 
Johnson-Cousins system, but comparison stars were chosen to reduce transformation effects 
whenever possible. 

For the reduction of the near-infrared (NIR) frames ESO's eclipse package was used 
(Devillard 2002). Photometric calibration of the ISAAC field was performed by the observa- 
tion of the UKIRT standard star FS 136 (Hawarden et al. 2001) with the VLT. Airmass cor- 
rections were applied according to the coefficients provided on ESO's Web pages (in units of 
mag airmass -1 ): kj=0.06, kn=0.06, A^=0.07. An independent calibration of the UKIRT K- 
band data, made using the standard star FS 130 (Hawarden et al. 2001), is in good agreement 
with the one described above. 3 Next we deduced the following magnitudes for the brightest 



2 IRAF is distributed by the National Optical Astronomical Observatories, which is operated by the 
Associated Universities for Research in Astronomy, Inc., under contract to the National Science Foundation. 

3 Following the detailed photometric work performed by Labbe et al. (2003) on the Hubble Deep Field 
South, we assumed that the ISAAC J s , H, K s filters match the faint NIR standard star system (Persson et 
al. 1998). We also refer the reader to Labbe et al. (2003; their section 2.1) for more details on the J s and 
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point-like object close to the optical transient (marked "Ref" in Fig. 1) at R.A. (J2000) = ll h 
33 m 04 s 27) Decl (J2000) = 25° 54' 22": J s = 17.10±0.10, H = 16.54±0.10, K s = 16.44±0.05 
4 . This calibration was confirmed with independent observations using the Calar Alto 3.5-m 
telescope equipped with the Omega Prime NIR camera on March 6, 2004. 

2.2. VLT spectroscopy 

The early report of the discovery of the GRB afterglow by Fox et al. (2003) allowed 
rapid spectroscopic observations of the optical transient, performed only 5 hrs after the burst, 
on 2003 Feb 26, 08:52-09:07 UT, with VLT Yepun/FORS2. A 900 s exposure using the 300V 
grism was obtained when the afterglow was at a magnitude of B = 19.8. Using a l'.'O slit, the 
3.3 A/pixel scale leads to a resolution of 13 A (FWHM, 870 km/s at A=4500 A) at O'.'9-l'.'l 
seeing. A second spectroscopic run was performed one night later, on 2003 Feb 27, 07:35- 
09:08 UT, using VLT Antu/FORSl equipped with the 600B grism. Five exposures, each 
one lasting 900 s, were taken at a time when the afterglow had a magnitude of B = 21.7. 
The 1.2 A/pixel scale gave a resolution of 7 A (FWHM, 470 km/s at A=4500 A) at l'.'O- 
L"2 seeing with a 1'.'3 slit. Flatfield and bias correction was applied in the usual way using 
IRAF. Background subtraction was performed and the wavelength calibration was done using 
HgCdHe and HgCdHe+Ar calibration lamps. The spectra were corrected for achromatic slit 
losses following the procedure described in Vreeswijk et al. (2004). The standard stars LTT 
6248 (spectral type A) and BD +33?2642 (spectral type B2 IV) were observed in order to 
flux-calibrate the spectra taken with FORS2 and FORS1, respectively. 

2.3. VLT polarimetry 

Optical polarimetric observations of the afterglow were performed with VLT Antu/FORSl 
and a Bessel R filter. Observations started on 2003 Feb 27.211 UT, approximately 1 day 
after the GRB, when the afterglow magnitude was about R = 20.7. Imaging polarimetry is 
achieved with a Wollaston prism splitting the image of each object in the field into the two 
orthogonal polarization components, which appear in adjacent areas of the CCD image. In 
this way, the measurement is insensitive of any variation in the source brightness during the 
observation. Two sets of images were acquired with four different angles (0°, 337?5, 315?0, 



K 8 filters. 

4 The optical magnitudes of this star are B = 20.15 ± 0.03, V = 19.07 ± 0.03, R c = 18.33 ± 0.13, I c = 
17.92 ±0.18 (Hcndcn 2003). 
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292?5) and with 15 min exposure each. Unfortunately, in one exposure a cosmic ray hit the 
afterglow, making any procedure unreliable. This set was therefore excluded from the anal- 
ysis. A polarimetric standard star, Hiltner652, was also observed in order to fix the offset 
between the polarization and the instrumental angles. Image reduction was performed using 
the eclipse tools (Devillard 2002), while photometry was extracted using the Gaia package 5 . 
The general procedure followed for FORS1 polarimetric observation analysis is extensively 
discussed, e.g., in Covino et al. (1999). 

2.4. Chandra X-ray observations 

A 40.17 ksec target-of-opportunity Chandra observation was initiated on 2003 February 
27 at 16:49 UT, about 1.54 days after the GRB. Only one chip was read out (the back- 
illuminated S3 chip) at a frame time of 0.44 sec. Three X-ray sources are detected in the 
8x8 arcmin 2 field of view with a significance of more than 5a, two of which are within 
the HETE/SXC error circle. The brightest source is detected with a mean count rate of 
0.010±0.001 cts s _1 , at R.A. = ll h 33 m 04:9; Decl. = +25° 53' 55" (J2000), coincident within 
the error (±0"5) with the position of the optical afterglow. This X-ray source also shows 
a fading by a factor of three during the observation, further strengthening its identification 
as the X-ray afterglow of GRB 030226. The other two sources are constant within their 
statistical errors, and detected at R.A. = ll h 33 m 02:8; Decl. = +25° 53' 55" and R.A. = 
U h 33 m 13 s 9 . Ded = +25 o 52 / 29", respectively. 

3. Results 

3.1. The optical/NIR/X-ray light curve 

In order to analyze the light curve of the GRB afterglow, we combined our data with pub- 
lished data from Gamma-Ray Burst Coordinated Network Circulars (GCNs 6 ; for references 
see Fig. 2) and Pandey et al. (2004). In doing so, we followed the standard procedure and 
fitted a broken power-law decay to all photometric data. Of the available multi-wavelength 
data, the i?-band light curve is sampled best. Therefore, we fitted the i?-band light curve, 
using the fitting equation suggested by Beuermann et al. (1999), but in the representation 



5 http : //star-www . dur . ac . uk/~pdraper/gaia/gaia . html 
6 http : //gen . gsf c . nasa . gov/gcn/gcn3 jir chive . html 
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given by Rhoads & Fruchter (2001): 

F„(t) = 2 l ' n F v {t b )[{t/t b )^ n + (t/t b y* n }- 1/n ■ (1) 

Here F v is the flux density 7 , a± and a 2 are the asymptotic decay indices before and after 
the light curve break, t is the time elapsed since the GRB trigger, t b the break time, and n 
the parameter which describes the smoothness of the break. In all cases we fitted apparent 
magnitudes after correction for Galactic extinction. The best fit of the combined i?-band 
data set is a x = 0.50 ± 0.35, a 2 = 2.55 ± 0.38, t b = 0.83 ± 0.10 days, and n = 0.91 ± 0.84 
(x 2 /d.o.f.= 2.8; Fig. 2). While similar values for a 2 have also been found for other afterglow 
light curves, the change in decay slope, a 2 — «i is rather large. If we use only our i?-band 
data (Table 2), the best fit is a x = 0.70 ± 0.25, a 2 = 2.66 ± 0.32, t b = 1.04 ± 0.12 days, and 
n — 1.2 ± 1.1 (x 2 /d.o.f.= 3.8), consistent with the previous results. Since no evidence for an 
underlying host was found in our data (§ 4.1), both fits assumed a magnitude for the host 
of m = 28 in all bands. In the following we will use the results of the first fit. 

Figure 4 shows the difference between observed and best fit magnitudes, adopting for 
all bands the light curve parameters obtained from the i?-band fit. Several small-scale 
fluctuations up to about 0.4 mag in both directions are apparent. In particular, on day 0.18 
the optical transient is notably brighter in K s than expected based on the light curve fit. 
On this day the afterglow was bright enough to do individual photometry of each NIR frame 
with consistent results compared to the stacked frames. This indicates that there were not 
any gross flattening errors since the optical transient was dithered across a fair fraction of 
the detector, so that some other explanation has to be invoked as to why this data point 
does not match the fit properly. Kulkarni et al. (2003) also noticed an unusual behavior of 
the afterglow light curve between 2.3 and 3.3 days after the burst. While flux fluctuations 
are also seen in the X-ray light curve observed with Chandra, a comparison of the strongest 
fluctuations in the R band (seen around 1.7 days after the GRB) with those in X-rays shows 
no obvious correlation between them (Fig. 5). However, in Fig. 5 four of the five i?-band 
data points were taken from the literature (GCNs), so that a potential small mismatch in 
the photometric calibration cannot be ruled out. 

3.2. The optical spectra 

The VLT spectra of the optical transient (Figs. 6, 7) show several absorption lines, 
but no prominent emission lines. The latter is not surprising, as the expected lines (e.g., 



7 We use the notation F v oc t a v 13 for the time and frequency dependence of the flux density. 
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[O II]A3727, [O III] A5007) are located outside our wavelength band, while Lya in emission 
is difficult to identify. The absorption lines are due to rest-frame ultraviolet metallic lines 
usually present in spectra of optical afterglows and high redshift galaxies (e.g., Castro et al. 
2003; Mirabal et al. 2003; Masetti et al. 2001, 2003; Savaglio et al. 2003; Shapley et al. 
2003). From these lines two redshift systems at z=1.962±0.001 and z=1.986±0.001 can be 
determined, consistent with earlier reports (Ando et al. 2003; Chornock & Filippenko 2003; 
Greiner et al. 2003a; Price et al. 2003b). In addition to the therein stated absorption line 
systems of Fell, Al II, and CIV, we find systems of OI, CII, Si IV, and Sill. In our first 
epoch spectrum we also detect a feature at 5711 and 5726A, which may be attributed to 
redshifted Mg II at ,2=1.042, a redshift system found via high-resolution spectroscopy with 
the Keck telescope (Chornock & Filippenko 2003; Price et al. 2003b). Based on our data 
we cannot identify the nature of this absorber (see also Fig. 11). Whether or not a similar 
feature at 5494 and 5517A as well as 6280 and 629lA is also due to redshifted Mg II is less 
certain. 

Following the detection of Fe II at a rest frame wavelength of 1608A (oscillator strength 
/=0.058; Prochaska et al. 2001) one would expect to see also Fell lines at rest frame 
wavelengths of 2344A (/=0.114), 2586A (/=0.069), 2600A (/=0.239) and perhaps even 
2374A (/=0.031). Although the lines are covered by the 300V spectrum of the first night, 
no clear line identification can be proposed as they blend with the telluric features at observer 
frame wavelengths around 7000 and 7800 A (Fig. 6). 

For the line fitting we used the IRAF/SPLOT task which assumes a Gaussian profile. 
Table 3 shows the line identifications together with the observer and rest-frame wavelengths, 
and the redshifts for each line detected. No line broadening and thus no velocity dispersion 
above the instrumental resolution is observed. In addition, the computed rest frame equiva- 
lent widths (EWs) are listed, which were calculated by dividing the measured observer frame 
EWs by (1 + z). Based on the measured equivalent width we calculated the column density, 
N, assuming the optically thin case (e.g., Spitzer 1968), 

where the symbols e,m e ,e , and c have their usual meaning for the SI system, and / is 
the oscillator strength of the corresponding transition. The resulting column densities are 
reported in Table 3. The majority of the observed absorption lines are strong (EW r >1 A) 
and may be saturated. In that case they deviate from the linear part of the curve of growth 
and the approximation given in Eq. (2) has to be considered as a lower limit for the derived 
column densities. 



Finally, the identified lines allow us to refine the redshift of the GRB to z—1. 986, or 
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larger. Assuming a flat universe with a matter density Qm = 0.3, vacuum density (cos- 
mological constant) = 0.7, and Hubble constant H = 65 km s -1 Mpc -1 , this redshift 
corresponds to a luminosity distance of 5.12 x 10 28 cm (16.6 Gpc) and a distance modulus 
of 46.1 mag. The look-back time is 11 Gyr. The difference in luminosity distance between 
both redshift systems is 245 Mpc, and the difference in velocity is 2410 km s _1 . 

3.3. The spectral energy distribution of the afterglow 

Of particular interest is the question whether or not the afterglow exhibited color evo- 
lution, as it has been reported for some GRBs (e.g., Bersier et al. 2003). However, given 
the relatively sparse data for GRB 030226 based on observing campaigns with different tele- 
scopes at different sites, we can only check for an obvious long-term trend in the color of 
the afterglow. Such a trend is not seen in our data (Fig. 3). Within the observational errors 
the break in the light curve was achromatic, and so was a potential flux variation on day 4.1 
(Fig. 2). This supports the view that the break was due to a collimated explosion (Rhoads 
1999; Meszaros & Rees 1999). On the other hand, around day 0.2 there is an excess flux in 
the K band (Fig. 4). Since the multi-color data were not taken simultaneously, we cannot 
decide whether this is indeed a color variation or caused by a flux variation as it has been 
seen in other well-studied afterglows (e.g., Lipkin et al. 2004). 

Based on the above discussion we conclude that the observed spectral energy distribution 
(SED) of the afterglow is best described by assuming no color variations. The shape of the 
light curve is then the same in all photometric bands. Using the magnitude-flux conversion 
according to Bessel (1979) and Bessel & Brett (1988), we find that the observed mean 
spectral slope, f3, across the BVRIJHK bands, corrected for Galactic extinction, is /3 — 
0.70 ± 0.03 (Fig. 8), with no evidence for additional reddening by dust in the GRB host 
galaxy (Ay (host) =0 mag). This holds for Milky Way-like dust as well as for SMC-like dust, 
for any assumed ratio of total-to-selective extinction between i?y=3 and 5. This result differs 
from Pandey et al. (2004) who based their conclusions exlusively on optical data obtained 
on days 0.62 and 1.79 after the burst (f3 ps 1). Pandey et al. argue that from the theoretical 
point of view the afterglow data are best understood if the intrinsic spectral slope was in 
fact ~0.55, requiring that the afterglow light was slightly reddened by dust in the GRB host 
galaxy. This spectral slope is indeed not much different from the one we deduce based on 
our larger data set, including NIR data, but without the need for additional reddening in 
the host galaxy. It is clear that the determination of f3 and the amount of reddening in a 
GRB host is very sensitive to the extent and quality of the data utilized. 

Is there evidence for dust in the host based on the X-ray data? The Chandra observation 
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took place after the break in the light curve (for details, see Pedersen et al. 2004b). The 
unabsorbed X-ray flux in the 0.3-10 keV band (source frame band of 0.9-30 keV) is 7.3 x 10" 14 
erg cm -2 s -1 , corresponding to a mean (isotropic) luminosity of 2.4xl0 45 erg s" 1 (averaged 
over the time interval of 1.54-2.04 days after the GRB). A free power-law fit to the X-ray 
spectrum with just local foreground absorption yields an energy index across the 0.3 to 10 
keV band of /9x=l-04±0.20 and a hydrogen column density iVfi = (4±1) xlO 20 cm" 2 (Fig. 9). 
Since the absorbing column is above the Galactic foreground value, we performed a second 
fit with two different absorbing columns, namely one fixed at the Galactic value (Nu(Ga\) 
= 1.8 xlO 20 cm" 2 ; Schlegel et al. 1998) and the other one at a redshift of the host galaxy 
of z = 1.986. The resulting energy index then is /3x=1.08±0.15 and the hydrogen column 
density in the host A^n(host) = (3.2±1.5) xlO 21 cm" 2 , corresponding to Ay(host) = 2.4 mag 
(for a Galactic dust-to-gas ratio). Evidence for so much extinction in the host galaxy is 
definitly not seen in the optical/NIR data, a discrepancy also known from other afterglows 
(Galama & Wijers 2001; Hjorth et al. 2003b; Stratta et al. 2004). 

Unfortunately, our VLT spectra cannot be used to derive the hydrogen column density 
in an independent way. The unfortunate position of Lja at the short edge of our VLT 
spectra and the uncertainty in the spectral shape longwards of the Lya-edge (Fig. 6) do not 
allow a quantitative treatment of the Lja absorption and therefore the hydrogen column 
density from the optical data. On the other hand, our polarization data also exclude a large 
amount of dust in the GRB host galaxy along the line of sight. 



3.4. The polarization of the afterglow 

Figure 10 shows the position of the afterglow in the Stokes Q-U- plane. The polarization 
level P is related to the Stokes parameters through P = \/Q 2 + U 2 , while the polarization 
angle is $ — \ arctan([//<3)- Our measurement yields for the afterglow Q = — 0.0028±0.0030, 
U = 0.0014 ± 0.0048 (la error), which translates, after correcting for the polarization bias 
(Wardle & Kronberg 1974), to a strict upper limit P < 1.1% (2a). This is one of the lowest 
limits ever reported for a GRB afterglow (e.g., Covino et al. 1999, 2003b; Rol et al. 2003; 
Barth et al. 2003), even if values as low as P ~ 0.5% have been measured for GRB 030329 
(Greiner et al. 2003b). The low polarization of the afterglow of GRB 030226 argues against 
substantial dust extinction in the GRB host along the line of sight, provided that the dust 
properties there are not much different from those of the dust in the local universe. We 
conclude that the polarimetry supports the view that the average dust-to-gas ratio in the 
GRB host galaxy along the line of sight is much lower than those in our Galaxy. 
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4. Discussion 

4.1. Limits on the host and a supernova component 

On our VLT images and in the afterglow light curve we do not find evidence for an 
underlying host galaxy (Fig. 11). Spectroscopic evidence for emission lines from an underly- 
ing host is missing too (although mainly due to the insufficient wavelength coverage of our 
spectra). Based on our latest i?-band imaging on March 13, i.e., 15 days after the burst, we 
conclude that i?host > 26.2. This estimate does not exclude the possibility that the host is 
indeed brighter, but some arcsec away from the burster (e.g., the objects 1 and 2 in Fig. 11), 
for which we have no spectroscopic confirmation. 

At the given redshift of z=1.986, we do not expect to see an underlying supernova (SN) 
component in our data: Assuming SN 1998bw as a template (for details, see Zeh et al. 
2004), this predicts an unabsorbed peak magnitude of R ~ 28 approximately at day 33 and 
a negligible contribution to the afterglow light within the first week after the burst. That 
neither an underlying host galaxy nor a SN component contributes substantially to the light 
of the optical transient, even 1 week after the burst, means that we essentially observed an 
uncontaminated GRB afterglow. 

4.2. The nature of the afterglow 

An achromatic break in the light curve is generally considered as evidence for a jetted 
explosion (Rhoads 1999; Meszaros & Rees 1999). Using standard procedures, for the given 
break time and redshift as well as gamma-ray fluence the beaming corrected energy output of 
GRB 030226 is consistent with the ensemble statistic of bursts; the results we obtain based 
on our deduced break time is basically the same as in Pandey et al. (2004), and we refer the 
reader to this paper for a more detailed discussion on this subject. 

A non-spherical explosion in the form of a highly relativistic jet should produce net 
polarization, and this is indeed what has been observed in a number of afterglows so far (for 
a review, see Covino et al. 2004). Geometrically simple models (e.g., Ghisellini & Lazzati 
1999; Sari 1999) assume a uniform jet, with a constant energy per unit solid angle accross 
the full jet. In this case, the polarization evolution has two peaks, separated by an epoch of 
null polarization, which roughly coincides with the break in the light curve. The polarization 
angle is expected to flip by 90° between the two peaks. More complex models (e.g., Rossi et 
al. 2004; Lloyd-Ronning et al. 2004) assume instead an energy distribution per solid angle 
which decreases towards the edge of the jet. In this case, even if the light curves are not 
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very different from the uniform case, the polarization is expected to show a peak close to 
the time of the break. Our polarization measurement of the afterglow of GRB 030226 was 
indeed performed very close to the break in the light curve. Interpreting this break as due 
to the jet effect, our low polarization value would thus favor a uniform jet over a structured 
jet. However, no solid conclusion can be drawn from a single upper limit, lacking monitoring 
of the polarization level and angle evolution with time. 

An inspection of Fig. 4 shows that all photometric data concentrate around the break 
time, while the other part of the light curve is more poorly sampled. Depending on the 
extent to which we merged data taken from different groups at different optical telescopes, 
our results vary between strong and weak evidence for a short-term variability in the light 
curve. These potential short-term fluctuations, primarily at later times, are reflected in the 
relatively large la error bars of the deduced light curve parameters a\ and a 2 . In particular, 
on day 4 after the burst, the afterglow flux is clearly below the expectations based on the 
light curve fit, and this holds for all optical bands (Fig. 4). Unfortunately, this is the only 
epoch after the break time where we have a homogeneous multi-color data set. Alternatively, 
there could be an excess of flux at later times (t > 5 days), which cannot be attributed to an 
underlying host or supernova, however. In this context it is not surprising that Pandey et 
al. (2004) derived a different post-break decay slope of a 2 = 2.05 ± 0.04 given the potential 
flux variations of the afterglow in combination with their much smaller monitoring time. 

A key parameter of the afterglow SED is the location of the cooling frequency u c , which 
separates fast-cooling from slow-cooling electrons (Sari et al. 1998). The combination of the 
Chandra X-ray data with the optical/NIR data places v c between the X-ray and the optical 
bands, in agreement with the conclusion drawn by Pandey et al. (2004). Thereby the spectral 
slopes are similar to those found for other afterglows (e.g., GRB 010222: Masetti et al. 2001; 
GRB 011211: Jakobsson et al. 2003). Standard fireball models for jetted explosions with 
p > 2 (Sari et al. 1998, 1999; Livio & Waxman 2000; Chevalier & Li 2000; Dai & Cheng 
2001) then predict p = 2(3 + 1 and a 2 = p (i.e., a 2 = 2.40 ± 0.06 for (3 = 0.70 ± 0.03; § 3.3), 
consistent with the observed late-time decay slope of the afterglow. 

Based on the above discussion we conclude that the afterglow of GRB 030226 showed 
all the properties already known from other bursts: an achromatic break around 1 day after 
the burst very likely due to a jetted outflow, a rapid decay thereafter, and possible short- 
term fluctuations. Finally, we can use the BOOTES-1 observations to constrain the optical 
emission before, during, and shortly after the genuine GRB event. A visual inspection of 
all BOOTES-1 frames reveals no such emission at the position of the optical afterglow, in 
particular simultaneously to the burst itself (Fig. 12). Therefore we derive an upper limit of 
R — 11.5 (due to cirrus clouds present in the sky at the time of the event) for the reverse 
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shock emission (if any) arising from this GRB, in contrast to the 9-th magnitude prompt 
flash emission recorded for GRB 990123 (Akerloff et al. 1999) which demonstrates that a 
bright, prompt optical emission is not a generic characteristic of GRBs. 

4.3. Clues on the GRB progenitor 

Of special interest are the highly ionized absorbers Si IV and C IV and the associated 
absorption line systems seen in our VLT spectra separated in velocity by 2400 km s _1 . 
Schaefer et al. (2003) and Mirabal et al. (2003) discussed a similar finding for the afterglow 
of GRB 021004 (z=2.3). The afterglow of this burst not only showed these lines, but four 
absorption line systems separated by 450, 990 and 3155 km s -1 (Mirabal et al. 2003; Castro- 
Tirado et al. 2004). Similar features were found in the afterglow of GRB 020813 (Barth et 
al. 2003). Schaefer et al. (2003) and Mirabal et al. (2003) have provided strong arguments 
that these lines, together with the observed absorption line systems, are likely to come 
from expanding shells around a massive Wolf-Rayet star (cf. Nugis & Lamers 2000; for a 
discussion, see also Wijers 2001; Chevalier et al. 2004). This could also explain the strong 
flux variations seen in some afterglows (e.g., GRB 000301C: Garnavich et al. 2000; Masetti 
et al. 2000; GRB 011211: Jakobsson et al. 2003; GRB 021004: Bersier et al. 2003; Holland 
et al. 2003; Lazzati et al. 2002; GRB 030329: Matheson et al. 2003), which might also be 
apparent for GRB 030226, even though with lower amplitude (Fig. 4). This is consistent 
with theoretical models that link long-duration GRBs to the explosions of massive stars (e.g., 
Fryer et al. 1999; Paczyhski 1998; Heger et al. 2003). 

For GRB 030226 the z— 1.986 system is clearly not a scaled version of the z— 1.962 
system. Basically, we identify two groups of ions. Those of Si IV, Fe II, and Al II have 
roughly the same rest-frame equivalent widths in both redshift systems, whereas those of 
the other species are approximately twice as abundant in the z= 1.986 system compared 
to the 2=1.962 system. Within the context of the Wolf-Rayet star model this would reflect 
different chemical compositions of different shells in the wind, a phenomenon that is observed 
for example in the outer ejecta of Eta Carina (Smith & Morse 2004). While evidence for 
'line locking' has been found in spectra taken from the afterglow of GRB 021004 (Savaglio et 
al. 2002; M0ller et al. 2002; see also the discussion in Mirabal et al. 2003), this phenomenon 
is not apparent in our spectra of the afterglow of GRB 030226. If not hidden because of the 
medium spectral resolution, this implies that the observed features in both absorption line 
systems are not due to radiative acceleration of the absorbing clouds. 

In principle, the presence of a strong Lja absorber in the spectrum of the afterglow 
of GRB 030226 implies a dense interstellar medium, possibly at the redshift of the burster. 
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However, in our spectra the redshift of the Ly a line cannot be accurately determined. We can 
only place the constraint 1.93< z <1.99. Similarly, a potential variability of the absorption 
lines is of interest since it could place the absorbing clouds close to the burster (e.g., Perna 
& Loeb 1998; Draine 2000). However, a comparison of our medium-resolution VLT spectra 
taken 0.2 and 1.2 days after the burst does not reveal evidence for a general variation in the 
line strengths. Unfortunately, the Lya line is difficult to analyze in this respect given its 
unfavorable position at the very blue side in our VLT spectra. 

While our spectral data support the view that the GRB progenitor was a massive star, 
there is no evidence for a stellar wind profile in the light curve data, similar to other cases 
(cf. Chevalier & Li 2000; Panaitescu & Kumar 2002). Within the framework of standard 
afterglow models with p > 2 the light curve parameter a.\ is related to the density profile of 
the circumburst medium, n(r) ~ r~ k , by oci — 3/4 (p — 1) + kj (8 — 2k) (e.g., Meszaros, Rees, 
& Wijers 1998; Dai & Wu 2003). In the ideal case k — means an ISM-shaped interstellar 
medium, while k = 2 refers to a stellar wind. Using a± = 0.50 ± 0.35 and p > 2 a density 
profile with k—2 is basically ruled out. Also the rather steep change in the decay slope 
between 0.5 and 2 days after the burst from a± = 0.50 ± 0.35 to «2 = 2.55 ± 0.38 might 
be difficult to explain with a standard wind profile. But does this really argue against WR 
stars as GRB progenitors? There is increasing observational evidence that the outer winds of 
these stars might indeed be very clumpy (e.g., Crowther et al. 2002) with the wind properties 
depending on the viewing direction, comparable to the one of Eta Carina (see Fig. 1 in van 
Boekel et al. 2003; Smith et al. 2003). An ansatz with k = 2 to describe the stellar wind 
might then be an oversimplification; in reality k might be notably smaller (Wijers 2001). 
More insight on the properties of winds from WR stars and more spectroscopic afterglow 
data are required in order to tackle this problem. A monitoring of the evolution of the cooling 
frequency across the optical/NIR bands in GRB afterglows, which might become possible in 
the Swift era (Gehrels 2004), might also help to clarify this issue. 



5. Summary 

We have presented a comprehensive data set of the afterglow of GRB 030226, consisting 
of optical, NIR, and X-ray observations, including VLT spectroscopy and VLT polarimetry. 
Concerning the nature of the GRB afterglow our basic findings are the following: 

I. The evolution of the afterglow was achromatic within the observational errors over the 
optical/NIR bands analogous to other bursts (e.g., GRB 000301C: Jensen et al. 2001; 
GRB 020405: Masetti et al. 2003). A break in the light curve occurred - 1 day after 
the burst. The break was rather smooth, similar to that found for, e.g., GRB 990510 
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(Harrison et al. 1999). 

2. Because the break was achromatic we interpret it as due to a jetted explosion. Close to 
the break time the intrinsic linear polarization of the afterglow was very low (< 1.1%), 
which favors a uniform jet over a structured jet. 

3. After correction for Galactic extinction, the best fit multi-color lightcurve has a spectral 
slope across the optical/NIR bands of /3 — 0.70 ± 0.03, with no sign of additional 
reddening due to dust in the host galaxy. This is in contrast to the Chandra X-ray 
data, which point to a relatively large hydrogen column density along the line of sight. 
This perhaps indicates a reduced dust-to-gas ratio in the GRB environment, similar to 
suggestions made for several other GRB afterglows (Galama & Wijers 2001; Hjorth et 
al. 2003b; Stratta et al. 2004). 

4. In the 0.3-10 keV region the spectral slope was (3 = 1.04 ± 0.20 on day 1.6-2, which 
together with the spectral slope across the optical/NIR bands places the cooling fre- 
quency between the optical bands and the X-ray region. 

In addition, based on BOOTES- 1 observations no prompt optical flash coincident with the 
gamma-ray burst was detected. We place an upper limit of R = 11.5 on the reverse shock 
emission arising from this burst. 

While our multi-wavelength observations trace the evolution of the afterglow, our VLT 
spectroscopy provides insight on the nature of the GRB environment, even though an under- 
lying host galaxy remained undetected down to a faint flux level (R > 26.2). In particular, 
beside a foreground absorber at 2=1.042 we found two absorption line systems at redshifts 
2=1.962 and 2=1.986, with the higher value as a lower limit for the redshift of the burster. 
Both line systems reveal an interstellar gas composed of highly ionized species (Si IV, CIV). 
This, together with the inferred separation in velocity of 2400 km s -1 between both systems, 
is very similar to observations of the afterglow of GRB 021004 ( Castro- Tirado et al. 2004; 
Schaefer et al. 2003; Mirabal et al. 2003), and possibly indicative of a fast stellar wind placed 
in a strong radiation field. Within this context the data support an association of GRBs with 
Wolf-Rayet stars. They also call for early spectroscopy of as many GRB afterglows as pos- 
sible in order to search for potential line features from stellar winds. It is surely premature 
to claim that all GRB progenitors are associated with WR stars, but among all progenitor 
models discussed so far in the literature (cf. Fryer et al. 1999) they are most promising from 
a theoretical point of view, and significantly supported by present-day observations. 
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A. AB offsets 

Table 1: Conversion constants for the transformation 
into the AB magnitude system (Oke & Gunn 1983). a 

Instrument Filter eff. wavelength AB off 



VLT/FORS1 


B Bess +34 


429 


nm 


-0.10 


VLT/FORS1 


\/_Bess +35 


554 


nm 


0.03 


VLT/FORS1 


i?_Bess +36 


657 


nm 


0.23 


VLT/FORS1 


J_Bess +37 


768 


nm 


0.44 


VLT/FORS2 


\/_Bess +75 


554 


nm 


0.06 


VLT/ISAAC 


Js 


1240 


nm 


0.90 


VLT/ISAAC 


H 


1650 


nm 


1.38 


VLT/ISAAC 


K s 


2160 


nm 


1.86 


NTT/EMMI 


B #605 


413.9 


nm 


-0.16 


NTT/EMMI 


V #606 


542.6 


nm 


0.02 


NTT/EMMI 


R #608 


641.0 


nm 


0.22 


NTT/SuSI2 


Bessel B #811 


421.2 


nm 


-0.06 


NTT/SuSI2 


Bessel R #813 


641.6 


nm 


0.21 


NOT/ALFOSC 


C/#7 


362 


nm 


0.72 


NOT/ALFOSC 


V #75 


530 


nm 


0.01 


TNG/Dolores 


C/#01 


361 


nm 


0.78 



a For the VLT/ISAAC filters the AB offset [mag] was taken from Labbe et al. (2003) and Gorosabel et al. 
(2003). For the other cases AB offsets were calculated in the present work using the effective wavelengths 
and detector quantum efficiencies as provided at the webpages of the corresponding observing sites. 
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Fig. 1.— VLT/ISAAC i^-band images of the afterglow of GRB 030226 obtained 4.5 hrs 
after the burst. The image size is approximately 100" x 100"; 'Ref denotes the reference 
star, which we used for the NIR photometry (see § 2.1). 
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Fig. 2. — Multi-color light curve of the optical transient following GRB 030226. Data 
from Table 2 are shown here (stars) together with data (open circles) reported in various 
Gamma-Ray Burst Circulars (Ando et al. 2003; Garnavich et al. 2003; Price & Warren 
2003; Guarnieri et al. 2003; von Braun et al. 2003; Rumyantsev et al. 2003a,b; Covino et 
al. 2003a; Fatkhullin et al. 2003; Semkov 2003) and by Pandey et al. (2004). Note that the 
i?-band data reported in GCNs # 1881, 1882, 1884, 1890, and 1908 were corrected for a 0.4 
mag offset at the zero point. 
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Fig. 3. — The color of the afterglow as a function of time, not corrected for Galactic extinc- 
tion. The color was calculated from AB magnitudes, with the AB offsets given in Table 1. 
Only data from telescopes with known AB offsets were used. Straight lines give the mean 
of the corresponding color. Even though we cannot resolve the origin for the offset between 
the NTT/EMMI and VLT/FORS2 V-band data during the first observing epoch, we favor 
the view that our data (Table 2) point to an achromatic evolution of the afterglow between 
0.2 and 4 days after the burst. 
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Fig. 4. — Observed minus best fit magnitudes (based on the .R-band fit; Fig. 2) for the data 
from Table 2 together with data reported in several GCNs (for references, see Fig. 2) and 
by Pandey et al. (2004). The scale at the left ordinate refers to U,V,I,H, the right to 
B, R, J, K . Am < means an excess of flux with respect to the fit. Symbols follow Fig. 2. 
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Fig. 5. — X-ray light curve obtained with Chandra (top panel) and the contemporaneous 
optical observations (bottom panel) from Fatkhullin et al. (2003), Semkov (2003), Pandey 
et al. (2004), and our observing run (Table 2). 
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Wavelength (A) 

Fig. 6. — Short-wavelength part of the VLT/FORS2 spectrum using the 300V grism taken on 
2003 February 26, five hours after the burst. Two distinct redshift systems at z— 1.962 and 
z=1.986 are detected. Line identifications are indicated. Note the broad Ly a absorption 
feature. Some telluric lines between ~6500 and 8000 A could not be removed during the 
calibration process. Also, the rather unexpected curvature in the blue wavelength range is 
probably unreal and may be a residual of the flux calibration. 
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Fig. 7. — VLT/FORS1 spectrum from the GRB afterglow using the 600B grism taken on 
2003 February 27, twenty nine hours after the burst and thus shortly after the break in the 
light curve. 
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Fig. 8. — Spectral energy distribution of the afterglow at the time of the light curve break 
(tb, § 3.1). The data are corrected for Galactic extinction. 
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Fig. 9. — The X-ray count spectrum (top) of the GRB 030226 afterglow, integrated over the 
full observation time (1.54-2.04 days after the GRB). A power-law model (full line) gives a 
reasonably accurate fit (x 2 /d.o.f.=0.7; deviations from the power law fit in units of a are 
shown in the bottom panel) with a spectral slope f3x = 1-04 ± 0.20. 
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Fig. 10. — Stokes Q-tZ-plane showing the position of the afterglow (big dot) compared with 
that of field stars. The average polarization of field stars is also close to zero, as it is expected 
towards such a low-reddening line of sight. 
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Fig. 11. — VLT/F0RS1 i?-band image taken on March 2, four days after the burst. At 
that time the magnitude of the optical afterglow was R—2A.2. There are two sources within 
2 / .'5 radius around the optical transient (indicated by numbers and arrows). Without having 
spectral data at hand, we cannot exclude that one of them is the host. However, in such a case 
the GRB would have occurred at a rather large projected distance from the luminosity center 
of the host: At 2=1.986 an angular distance of 1" corresponds to 9 kpc. It is conceivable 
that one of these sources is responsible for the foreground absorber of Mg II at z— 1.042 we 
found in our spectra (§ 3.2). 
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Fig. 12. — Unfiltered exposures (180-s each) obtained by the BOOTES-1 wide field camera 
covering the GRB 030226 error box. The images covered the period 03:41 - 03:52 UT on 26 
Feb 2003 (i.e., between 5.3 min before the trigger and 5.8 min after the trigger). The ticks 
mark the position of the optical afterglow. The despicted field is 2x2 deg 2 . Note that there 
is some degree of geometrical distortion due to the fact that the images were taken with a 
wide- field optical system. North is up and east to the left. 
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Table 2 — Continued 
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Table 2 — Continued 
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Table 3: Line identifications 3. 
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"Symbols: oscillator strengths / (Prochaska et al. 2001), observer-frame wavelengths A obs , corresponding 
absorption redshifts z a ^ Sl rest-frame equivalent widths (EW r ) and deduced column densities (iV[cm~ 2 ]) of 
the absorption lines detected in VLT spectra. Typical errors in redshift are ±0.001. Column densities were 
calculated for the optically thin case. Note that the C II and C IV doublet is not resolved in our spectra. 
The column densities for these ions refer to the transition with the largest oscillator strength and neglect 
the other component. 



